Measurements of Kondo and spin splitting in single-electron transistors 
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We measure the spin splitting in a magnetic field B of localized states in single-electron transistors 
using a new method, inelastic spin-flip cotunneling. Because it involves only internal excitations, 
this technique gives the most precise value of the Zeeman energy A — IqIubB. In the same devices 
we also measure the splitting with B of the Kondo peak in differential conductance. The Kondo 
splitting appears only above a threshold field as predicted by theory. However, the magnitude of 
the Kondo splitting at high fields exceeds 2\g\iJ,BB in disagreement with theory. 
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The Kondo state, formed when conduction electrons 
screen a magnetic impurity, has received new attention 
since its observation in nanostructures, because it is a 
spin-entangled state of the many-electron system. While 
the equilibrium properties of Kondo systems are de- 
scribed quantitatively by theory, non-equilibrium Kondo 
phenomena have proved harder to understand. Single- 
electron transistors (SETs) provide the unique possibility 
of exploring these non-equilibrium phenomena. An SET 
consists of a confined droplet of electrons, called an arti- 
ficial atom or a quantum dot, coupled by tunnel barriers 
to two conducting leads, called the source and the drain. 
The electrochemical potential of the dot, as well as the 
coupling between the dot and the leads, can be tuned 
by changingvoltages on electrodes. The Kondo effect in 
SETs 111 occurs because a dot in a non-zero spin 

state coupled to its leads is analogous to a magnetic im- 
purity coupled to the electrons in a host metal. Kondo 
correlations develop because the dot spin is screened by 
the spins of electrons in the leads. In an SET the Kondo 
state can be studied out of equilibrium by applying a DC 
voltage Vds between the source and the drain. 

Non-equilibrium Kondo physics is probed by the mea- 
surement of the splitting of the Kondo peak in differential 
conductance as a function of magnetic field B 0, 0, 0, ■ 
Meir et al. Q predict that in a magnetic field the Kondo 
peaks occur at eVds = ±A. Here A = \g\iji,BB is the 
Zeeman energy for spin splitting and /ub = 58 /ieV/T 
is the Bohr magneton. Cronenwett et al. have mea- 
sured the splitting of a Kondo peak and have found good 
agreement with this prediction, using g = —0.44 for bulk 
GaAs. However, more recent calculations by Costi 6] 
predict that the Kondo peak splitting should appear only 
above a critical magnetic field, and Moore and Wen 7] 
predict that the splitting should be smaller than 2A at 
all fields. 

Here we report measurements comparing the Kondo 
peak splitting to A in SETs. Using electron addition 
spectroscopy and a new, more precise method, inelastic 
spin-flip cotunneling, we determine A and hence the g- 
factor of our SETs. We find that the Kondo splitting 
appears only above a critical field as predicted by theory. 



However, the magnitude of the Kondo splitting at high 
fields is noticeably larger than 2 A, contrary to theoretical 
predictions. 

The two SETs we have studied are similar to those used 
by Goldhaber-Gordon et al. and details about the 

AlGaAs/GaAs heterostructure can be found in the lat- 
ter references. A two dimensional electron gas (2DEG) 
is formed at the AlGaAs/GaAs interface with an elec- 
tron density of 8.1 x 10^^ cm~^ and a mobility of 10^ 
cm^/Vs at 4.2 K. Magneto-transport shows that only 
one sub-band of the 2DEG is occupied. Electron-beam 
lithography is used to define the gate electrodes shown 
in the inset of Fig. ^c). Applying a negative voltage 
to these electrodes depletes the 2DEG underneath them 
and forms an artificial atom of about 50 electrons iso- 
lated by two tunnel barriers from the remaining 2DEG 
regions, the source and drain leads. The voltage on the 
gate electrode g is denoted Vg. 

The two SETs were studied in different dilution refrig- 
erators. SET 1 was measured in a Leiden Cryogenics 400 
/iW dilution refrigerator with a 14 T magnet and a low- 
est electron temperature of 20 — 30 mK. The sample was 
aligned so that the 2DEG was parallel to the magnetic 
field to better than 0.5 degree. SET 2 was measured in a 
75 /iW Oxford Instruments dilution refrigerator with an 
8 T magnet and a lowest electron temperature of about 
110 mK. This sample could only be aligned parallel to 
the magnetic field to within a few degrees. To measure 
the differential conductance dl/dVds , we added a small 
sinusoidal modulation to Vds and measured the resulting 
current with a preamplifier and a lock-in amplifier. 

One way to measure the Zeeman splitting of orbital 
levels in an SET is electron addition spectroscopy |lol . lll| . 
Figure nia) shows a plot of dl/dVds as a function of Vds 
and Vg while Fig. ^b) shows dl/dVds at the value of Vg 
marked on Fig. QJa) with a dashed line, for which the 
occupancy of the dot is even. As we increase Vds at this 
fixed Vg, we lower the Fermi energy of the drain lead. 
The first peak in dl/dVds occurs when the Fermi energy 
of the drain aligns with the higher energy state of the 
spin split orbital state in the SET. Since even occupancy 
implies both spin states are filled at Vds = 0, a second 
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FIG. 1: (a) Differential conductance as a function of Vds and 
Vg at B = 7.25 T for SET 2. The odd and even Coulomb 
valleys are labeled, (b) Differential conductance as a function 
of drain-source voltage at the gate voltage marked by the 
dashed line in (a). The dots denote data while the solid line 
shows a fit to the sum of two cosh~^ functions, (c) Zeeman 
splitting as a function of magnetic field determined by fitting 
curves like that shown in (b) . The solid line shows a fit to the 
data to determine \g\. The inset shows an electron micrograph 
of an SET similar to those we studied. 



peak appears when the lower energy spin state enters 
the transport window. The separation of the two peaks 
is A/e(l — ad), where ad is the ratio of the dot-drain 
capacitance to the total capacitance. The (1 — ad) factor 
arises because applying a positive voltage to the drain 
slightly lowers the energies of the levels in the dot. 

To find the peak splitting we fit data like those in 
Fig. nib) to the sum of two cosh~^ functions; the spe- 
cific form chosen does not affect the measurement of the 
splitting. Using the slopes of a Coulomb-blockade dia- 
mond, part of which is shown in Fig. Ql^a) , we determine 
ad = 0.15 ± 0.02 and hence A. Figure ^c) shows A as 
a function of _B. At each field we fit several traces like 
that in Fig. Q^b). The scatter in the data arises because 
charge fluctuations near the SET are equivalent to gate 
voltage fluctuations and cause fluctuations in the peak 
positions. Fitting the data in Fig. ^c) to a line gives 
|5| =0.18 ±0.04. 

We have developed a new, more precise method of mea- 
suring A using inelastic spin-flip cotunneling. Although 
resonant tunneling is prohibited in the Coulomb blockade 
regime, an electron can tunnel through the SET via a vir- 
tual intermediate state, a process known as cotunneling 



[l^. Cotunneling can proceed elastically and leave the 
SET in its ground state, or inelastically, leaving it in an 
excited state. De Franceschi et al. observed inelastic 
cotunneling through an excited orbital state of an SET 
with an even number of electrons. In our experiments the 
ground state has an odd number of electrons and the role 
of the excited state is played by the higher energy spin 
state. The elastic and inelastic cotunneling processes are 
illustrated in Fig.|2a)-(d). If |eVds| < A only the elastic 
process is possible, but if |eVds| > A both the elastic 
and inelastic processes are possible and we expect steps 
in dl/dVds at |eVds| = A. 

This behavior is illustrated in Fig.l^fa) and (b). Note 
that the spin-flip cotunneling gap is independent of gate 
voltage and, as shown in Fig. ^a), is only observed in 
a magnetic field. This is in contrast to cotunneling fea- 
tures associated with excited orbital states, like those in 
the Coulomb valley on the far left in Fig. 3(a), which 
depend on Vg because changing the gate voltage alters 
the shape of the confining potential 0, 0| . The orbital 
cotunneling features are also present at B = 0. To our 
knowledge, this is the first observation of the threshold 
for spin- flip cotunneling in an SET. 

Increasing the temperature broadens the inelastic co- 
tunneling threshold. Lambe and Jaklevic IS] have shown 
that for an inelastic process with negligible intrinsic 
width the lineshape is given by 
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In this equation, Ag is the conductance from elastic co- 
tunneling and Ai describes the additional contribution 
from inelastic cotunneling. is a function defined by 
F{x) = [1 + (a: — 1) exp(a;)]/[exp(x) — 1]^. This hneshape 
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FIG. 2: Cotunneling processes through a Zeeman split orbital, 
which can occur when the dot is in the Coulomb blockade 
regime with an odd number of electrons as in (a). We take 
spin up to be the lower energy state. A spin-down electron 
from the lead can tunnel onto the dot forming the virtual 
intermediate state shown in (b). For any value of Vds, the 
spin-down electron can tunnel off the dot as shown in (c), 
resulting in elastic cotunneling. If |eVds| > A = \g\fiBB, 
the spin-up electron can also tunnel off the dot resulting in 
inelastic cotunneling as shown in (d). 
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FIG. 3: (a) Differential conductance as a function of Vds and 
Vg at B = 8 T for SET 1. The inelastic spin-flip cotunneling 
gap is in the middle of the center diamond and is marked by 
arrows, (b) Differential conductance as a function of drain- 
source voltage in the middle of the Coulomb valley for the 
same sample as in (a), but with different voltages on the elec- 
trodes and B = 11 T. The dots are the data and the solid 
hne is a fit described in the text, (c) Temperatures extracted 
from fitting the cotunneling gap of SET 2 at B = 7.65 T. The 
horizontal axis is the mixing chamber temperature Tmc and 
the vertical axis is the temperature Tpit extracted from fitting 
the gap as discussed in the text. The line shows Tmc = Tpit. 



has steps centered at Vds = ±A/e with width 5.4/cBT/e. 
A fit to Eq. ^is shown as the solid line in Fig.lSJb) and is 
in excellent agreement with the data. Figure |5Jc) shows 
the temperature Tpit extracted from fits to data taken at 
different temperatures with SET 2. The sharpest cotun- 
neling step that we have observed (shown in Fig. Olb)) 
is for SET 1 and has width 22 fieV. This can be used to 
place a lower bound on the spin decoherence time because 
the intrinsic width of the cotunneling step is inversely 
proportional to the decoherence time. 

Figure Efa) shows the evolution with field of a cotun- 
neling feature in SET 2. The splittings A from fits to Eq. 
^are plotted as functions of B in Fig.^Jb). Fitting these 
data to a line gives \g\ = 0.16 ± 0.02, consistent with the 
value from addition spectroscopy but with greater preci- 
sion. Measurements for SET 1 in high magnetic fields, 
shown in Figs. Efc) and (d), give \g\ = 0.152 ± 0.006, 
consistent with the value for SET 2. The inelastic cotun- 
neling method, which measures the spectrum of internal 
excitations of the dot, is more precise than electron addi- 
tion spectroscopy for two reasons. First, the cotunneling 
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FIG. 4: (a) Evolution of a cotunneling gap (dots) and a Kondo 
peak (solid line) with increasing magnetic field in SET 2. The 
vertical scale on the left (right) is for the cotunneling data 
(Kondo data). For the cotunneling data, the coupling to the 
leads was reduced so that no Kondo peak was seen at zero 
field at the measurement temperature T ~ 110 mK. The solid 
lines through the cotunneling data at high fields are examples 
of fits to equation The dashed lines at the three highest 
fields mark the center of the cotunneling edge; the lines are 
2 A apart, (b) Splittings from cotunneling gaps (open circles) 
and Kondo peaks (solid squares) from SET 2. The solid line 
through the cotunneling data is a linear fit. The line through 
the Kondo data is the result of fixing the slope at the value 
obtained from the cotunneling data and fitting the data above 
3 T to extract the intercept, (c) Evolution of a Kondo peak 
into a cotunneling gap with increasing magnetic held in SET 
1. The solid line in the 5 T trace shows a fit as described 
in the text, (d) Splittings from the cotunneling gaps and 
Kondo peaks for SET 1. The Kondo splitting is determined 
at low field and the cotunneling gap at high field, but all 
other parameters are unchanged. We fit the data with lines 
the same way as for (b). 



spectrum is independent of the chemical potential of the 
dot, making the splitting insensitive to small charge fluc- 
tuations near the SET. Second, the cotunneling method 
is inherently more precise because its intrinsic width is 
smaller than the width of the charging peak, which is 
determined by the coupling to the leads and is approxi- 
mately r = aS^ieV for SET 2 in Fig. [T] 

The splitting of the Kondo peak at high fields is larger 
than the cotunneling gap, as can be seen, for example, 
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from the data in Fig. QJa). To compare our findings to 
theory, we foUow Moore and Wen |7| who identify the 
field-induced sphtting between the two peaks in their 
calculated spectral function with the splitting in Vds be- 
tween the peaks in dl/dVds- We define Ak/c as half 
the separation in Vds between the two peaks. Figure 
2fb) shows Ak for SET 2 from data like those shown in 
Fig. EJa). There is a clear offset between Ak and the 
cotunneling data at high fields; fitting the data above 
3 T to a line with \g\ fixed at the value obtained from 
cotunneling gives an offset of f 2 ± 1 /iteV. 

The Kondo peak evolves smoothly into a cotunneling 
threshold with increasing B when there is an odd number 
of electrons in the dot. Figure 0fc) shows this evolution 
for SET f . Note that the values of A in Fig. Efd) are 
extracted from the high-field data for which the Kondo 
peak is absent. The Ak are extracted from the low- field 
data and are plotted in Fig.^d). We fit the data to a hue 
with the slope fixed at the value obtained from the cotun- 
neling data points, \g\ = 0.152. The fit gives a y-intercept 
of 10 ± 2 iieV, similar to the result for SET 2. Note that 
in contrast to the Kondo data, both the addition spec- 
troscopy and cotunneling measurements extrapolate to 
zero within the errors at B = 0. Contrary to theoretical 
predictions 0|, Ak is larger than \g\fiBB. 

Another feature of the Kondo data is that the splitting 
only appears above a threshold magnetic field. For ex- 
ample, the 1.25 T data in Fig. 2Ja) shows no evidence of 
splitting. The data for Ak in Fig. ^b) shows a thresh- 
old for splitting near 2 T, which is where \g\iJ,BB/kB ~ 
Tk ~ 300 mK. Costi predicts a threshold field for the 
Kondo splitting. In this respect our results are consistent 
with theory. 

The value of |g| in our SETs is 0.16, much smaller than 
the bulk GaAs value of 0.44. The small g-factor may arise 
because the electron density in our SETs is much larger 
than in typical devices. Variations in the g-factor with 
density have previously been observed in AlGaAs/GaAs 
heterostructures [ll[ll[ll[ll|23. We beheve our smaU 
value results from penetration of the electrons into the 
Alo.3Gao.7As which has g = -1-0.4, and the large 2DEG 
density in our SETs enhances this effect. Linearly extrap- 
olating the results of Jiang and Yablonovitch 'l9| gives a 
change in \g\ that is about half of that we observe. We 
plan a more realistic calculation to test this hypothesis 
more carefully. 

The observation that Ak > l^l/is-B is surprising. 
Moore and Wen predict that Ak < \g\fJ,BB at all 
B, whereas we find that Ak > \g\nBB by « 10 /icV at 
high fields. Our results differ from those of Goldhaber- 
Gordon et al. Ij. However, the measurements reported 
by Goldhaber-Gordon et al. are taken with the magnetic 
field normal to the 2DEG and could involve effects from 
the orbital motion of the electrons. Such effects are not 
considered by Moore and Wen and are eliminated in 
the present experiments by orienting the 2DEG parallel 



to the magnetic field. One might ask if elastic cotunnel- 
ing forms a background that could explain our observa- 
tions. Using zero field data like those in Fig. QJa), only 
over a wider range of Vds, we fit the cotunneling conduc- 
tance outside the Kondo peak region to a parabola and 
subtract this from our Kondo splitting data from SET 2. 
This decreases the high-field offset by about 5%, which 
does not affect our conclusions. However, such a sub- 
traction is not well motivated since elastic and inelastic 
cotunneling processes contribute coherently to the Kondo 
correlations. Thus the observation that Ak > jgl/isi? is 
not understood and suggests that more theoretical anal- 
ysis of the non-equilibrium Kondo problem is important. 
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